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STOCHASTIC EQUATIONS, FLOWS 
AND MEASURE- VALUED PROCESSES 

By Donald A. Dawson 1 and Zenghu Li 2 

$h ■ Carleton University and Beijing Normal University 

We first prove some general results on pathwise uniqueness, com- 
parison property and existence of nonnegative strong solutions of 
stochastic equations driven by white noises and Poisson random mea- 
sures. The results are then used to prove the strong existence of 
two classes of stochastic flows associated with coalescents with mul- 
tiple collisions, that is, generalized Fleming- Viot flows and flows of 
continuous-state branching processes with immigration. One of them 
unifies the different treatments of three kinds of flows in Bertoin and 
Le Gall [Ann. Inst. H. Poincare Probab. Statist. 41 (2005) 307-333]. 
Two scaling limit theorems for the generalized Fleming- Viot flows 
are proved, which lead to sub-critical branching immigration super- 
processes. From those theorems we derive easily a generalization of 
the limit theorem for finite point motions of the flows in Bertoin and 
Le Gall [Illinois J. Math. 50 (2006) 147-181]. 

1. Introduction. A class of stochastic flows of bridges were introduced 
by Bertoin and Le Gall (2003) to study the coalescent processes with mul- 
tiple collisions of Pitman (1999) [see also Sagitov (1999)]. The law of such 
a coalescent process is determined by a finite measure A(dz) on [0,1]. The 
Kingman coalescent corresponds to A = 5o and the Bolthausen-Sznitman co- 
alescent corresponds to A = Lebesgue measure on [0, 1] [see Bolthausen and 
Sznitman (1998) and Kingman (1982)]. In fact, Bertoin and Le Gall (2003) 
established a remarkable connection between the coalescents with multiple 
collisions and the stochastic flows of bridges. Based on this connection, they 
have developed a theory of the coalescents and the flows in the series of 
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papers [see Bertoin and Le Gall (2003, 2005, 2006)]. We refer the reader to 
Le Jan and Raimond (2004), Ma and Xiang (2001) and Xiang (2009) for the 
study of stochastic flows of mappings and measures in abstract settings. 

Let {B s j ■ —oo < s <t < oo} be the stochastic flow of bridges associated 
to a A-coalescent in the sense of Bertoin and Le Gall (2003). A number of 
precise characterizations of the flow {B_t,o{v) '-t > 0,v G [0,1]} were given 
in Bertoin and Le Gall (2003). For any t > 0, the function v t— > B_t,o( v ) 
induces a random probability measure pt(dv) on [0, 1]. The process {pt :t > 
0} was characterized in Bertoin and Le Gall (2003) as the unique solution of 
a martingale problem. In fact, this process is a measure-valued dual to the 
A-coalescent process. It was also pointed out in Bertoin and Le Gall (2003) 
that {pt :t > 0} can be regarded as a generalized Fleming-Viot process [see 
also Donnelly and Kurtz (1999a, 1999b)]. 

Let A(dz) be a finite measure on [0, 1] such that A({0}) = 0, and let 
{M(ds, dz, du)} be a Poisson random measure on (0, oo) x (0, l] 2 with inten- 
sity z~ 2 ds A(dz) du. It was proved in Bertoin and Le Gall (2005) that there 
is weak solution flow {Xt(v) : t > 0, v G [0, 1]} to the stochastic equation 

(1.1) X t (v) = v+ [ [ [ z[l {u < Xa _ {v)} -X s -(v)]M(ds,dz,du). 

Jo Jo Jo 

Moreover, Bertoin and Le Gall (2005) showed that for any < r% < ■ ■ ■ < 
r p < 1 the p-point motion {(B_ t ,o( r i), ■ ■ ■ , B^ t ,o( r p)) '■ t > 0} is equivalent 
to {(X t (ri), . . . ,X t (r p )) :t > 0}. Therefore, the solutions of (1.1) give a re- 
alization of the flow of bridges associated with the A-coalescent process. 
A separate treatment for the Kingman coalescent flow was also given in 
Bertoin and Le Gall (2005). In that case they showed the p-point motion 
{(-B_ t] o(ri), . . . ,B-t,o(fp)) - t>0} is a diffusion process in 

D p ■= {x = (xi, . . . , x p ) G W : < x x < ■ • • < x p < 1} 

with generator Aq defined by 

(1.2) Aof(x) = 2^2 W 1 " x ^j) dx . ( x )- 

i,j = l 1 3 

Given a A-coalescent flow {B s j : — oo < s < t < oo}, we define the flow of 
inverses by 

B;j(v)=mi{u£ [0,l]:B s>t (u)>v}, v G [0, 1), 

and B~l(l) = BJ$(1-). In the Kin gman coalescent case, it was proved in 

Bertoin and Le Gall (2005) that the p-point motion { {Bq\ (r\ ) , . . . , BqI(t p )) : 
t > 0} is a diffusion process in D p with generator A± given by 

(1.3) Axfix) = A f(x) + J2( \ ~ xi) !£(*), 
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where Aq is given by (1.2). The analogous characterization for the A-coa- 
lescent flow with A({0}) = was also provided in Bertoin and Le Gall (2005). 
Those results give deep insights into the structures of the stochastic flows 
associated with the A-coalescents. 

The asymptotic properties of A-coalescent flows were studied in Bertoin 
and Le Gall (2006). For each integer k > 1 let A k (dx) be a finite measure on 
[0,1] with A fe ({0}) = and let {X k (t, v) : t > 0, v G [0,1]} be defined by (1.1) 
from a Poisson random measure {Mf c (ds,dz,du)} on (0, oo) x (0, l] 2 with 
intensity z~ 2 ds A; c (dz) du. Suppose that z~ 2 (z A z 2 )A^.(A;~ 1 dz) converges 
weakly as k — > oo to a finite measure on (0, oo) denoted by z~ 2 (z A z 2 )A(dz). 
By a limit theorem of Bertoin and Le Gall (2006) the rescaled p-point mo- 
tion {(kX/ ( .(kt,ri/k), . . . ,kXk(kt,r p /k)) :t > 0} converges in distribution to 
that of the weak solution flow of the stochastic equation 



(1.4) Y t (v)=v+ / / xl {u < Ys „(v)}N(ds,dx,du), 

Jo Jo Jo 

where N(ds,dx,du) is a compensated Poisson random measure on [0,oo) x 
(0,oo) 2 with intensity z~ 2 ds A{dz) du. It was pointed out in Bertoin and 
Le Gall (2006) that the solution of (1.4) is a special critical continuous- state 
branching process (CB-process). 

In this paper we study two classes of stochastic flows defined by stochas- 
tic equations that generalize (1.1) and (1.4). We shall first treat the gen- 
eralization of (1.4) since it involves simpler structures. Suppose that a > 
and b are constants, v \— > ^y(v) is a nonnegative and nondecreasing continu- 
ous function on [0,oo) and (z A z 2 )m(dz) is a finite measures on (0,oo). Let 
{W(ds,du)} be a white noise on (0,oo) 2 based on the Lebesgue measure 
dsdu. Let {N (ds , dz , du)} be a Poisson random measure on (0,oo) 3 with 
intensity dsm(dz) du. Let {N (ds , dz , du)} be the compensated measure of 
{N (ds , dz , du)} . We shall see that for any v > there is a pathwise unique 
nonnegative solution of the stochastic equation 

ft pY a —(v) ft 

Y t (v)=v + o I / ' W(ds,du)+ / [j(v) -bY s _(v)]ds 

(1.5) J ° J ° J ° 

r-t roo rY 3 -(v) 

+ ij I zN(ds,dz,du). 



J0 J0 JO 

It is not hard to show each solution Y(v) = {Yt(v) :t > 0} is a continuous- 
state branching process with immigration (CBI-process). Then it is natural to 
call the two-parameter process {Yt(v) :t > 0,v > 0} a flow of CBI-processes. 
We prove that the flow has a version with the following properties: 

(i) for each v > 0, 1 1- > Y t (v) is a cadlag process on [0, oo) and solves (1.5); 

(ii) for each t > 0, «h> Yt(v) is a nonnegative and nondecreasing cadlag 
process on [0, oo). 
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The proof of those properties is based on the observation that {Y{v) : v > 0} 
is a path-valued process with independent increments. For any t > 0, the ran- 
dom function t> i— > Yt(v) induces a random Radon measure Yt(dv) on [0,oo). 
We shall see that {Yt '■ t > 0} is actually an immigration superprocess in the 
sense of Li (2011) with trivial underlying spatial motion. One could replace 
the diffusion term in (1.5) by the stochastic integral a j Q y/Y s -{v) dW(s) 
using a one-dimensional Brownian motion :t > 0} as in Dawson and 

Li (2006). The resulted equation defines an equivalent CBI-process for any 
fixed v > 0, but it does not give an equivalent flow. 

To describe our generalization of (1.1), let us assume that a > and b > 
are constants, v i— > 7(f) is a nondecreasing continuous function on [0, 1] such 
that < 7(u) < 1 for all < v < 1 and z 2 v(dz) is a finite measure on (0, 1]. 
Let {B(ds,du)} be a white noise on (0,oo) x (0,1] based on dsdu, and 
let {M(ds,dz,du)} be a Poisson random measure on (0,oo) x (0, l] 2 with 
intensity dsv(dz)du. We show that for any v G [0,1] there is a pathwise 
unique solution X{v) = {Xt(v) : t > 0} to the equation 



Clearly, the above equation unifies and generalizes the flows described 
by (1.1), (1-2) and (1.3). Here it is essential to use the white noise as the dif- 
fusion driving force. We show there is a version of the random field {Xt(v) : 
t > 0, < v < 1} with the following properties: 

(i) for each v £ [0, 1], tt-t X t (v) is cadlag on [0,oo) and solves (1.6); 

(ii) for each t > 0, v h-> Xt(«) is nondecreasing and cadlag on [0, 1] with 
X t (0) > and X t (l) < 1. 

We refer to {X^(t>) :t > 0,0 < v < 1} as a generalized Fleming-Viot flow 
following Bertoin and Le Gall (2003, 2005, 2006). In particular, our result 
gives the strong existence of the flows associated with the coalescents with 
multiple collisions. The study of this flow is more involved than the one 
defined by (1.5) as the path-valued process {X(v) : < v < 1} does not have 
independent increments. However, we shall see it is still an inhomogeneous 
Markov process. From the random field {Xf(v) :t > 0,0 < v < 1} we can 
define a cadlag sub-probability- valued process {Xt : t > 0} on [0,1], which is 
a counterpart of the generalized Fleming-Viot process of Bertoin and Le Gall 
(2003). We prove two scaling limit theorems for the generalized Fleming- 
Viot processes, which lead to a special form of the immigration superprocess 




(1.6) 
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defined from (1.5). Prom the theorems we derive easily a generalization of 
the limit theorem for the finite point motions in Bertoin and Le Gall (2006). 

The techniques of this paper are mainly based on the strong solutions 
of (1.5) and (1.6), which are different from those of Bertoin and Le Gall 
(2005, 2006). In Section 2 we give some general results for the pathwise 
uniqueness, comparison property and existence of nonnegative strong solu- 
tions of stochastic equations driven by white noises and Poisson random 
measures. Those extend the results in Fu and Li (2010) and provide the 
basis for the investigation of the strong solution flows of (1.5) and (1.6). 
They should also be of interest on their own right. In Section 3 we study 
the flows of CBI-processes and their associated immigration superprocesses. 
The generalized Fleming- Viot flows are discussed in Section 4. Finally, we 
prove the scaling limit theorems in Section 5. 

Notation. For a measure (i and a function / on a measurable space (E, S) 
write (fi, f) = f E f dfi if the integral exists. For any a > let M[0, a] be the 
set of finite measures on [0, a] endowed with the topology of weak conver- 
gence. Let Mi[0,a] be the subspace of M[0,a] consisting of sub-probability 
measures. Let -B[0,a] be the Banach space of bounded Borel functions on 
[0,a] endowed with the supremum norm || • ||, and let C[0,a] denote its sub- 
space of continuous functions. We use -B[0,a] + and C[0,a] + to denote the 
subclasses of nonnegative elements. Throughout this paper, we make the 
conventions 

["=[ and r=f 

J a J (a,b] J a J (a,oo) 

for any b > a > 0. Given a function / defined on a subset of M, we write 

A z f(x) = f(x + z)-f(x) and D z f(x) = AJ(x)-f'(x)z 

for x,z £M if the right-hand side is meaningful. Let A denote the Lebesgue 
measure on [0,oo). 

2. Strong solutions of stochastic equations. In this section, we prove 
some results on stochastic equations of one-dimensional processes driven 
by white noises and Poisson random measures. The results extend those of 
Fu and Li (2010). Since our aim is to apply the results to the generalized 
Fleming- Viot flows and the flows of CBI-processes, we only discuss equations 
of nonnegative processes. However, the arguments can be modified to deal 
with general one-dimensional equations. 

Let E, Uq and U\ be separable topological spaces whose topologies can 
be defined by complete metrics. Suppose that ir(dz), fJ-o(du) and Hi(du) are 
cr-finite Borel measures on E, Uq and U\ , respectively. We say the parame- 
ters (cr,b,go,gi) are admissible if: 

• x i — y b{x) is a continuous function on R + satisfying 6(0) > 0; 
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• (x,u) \— > o~(x,u) is a Borel function on M + x E satisfying a(0,u) = for 
u £ E; 

• (x, u) i — y go(x, u) is a Borel function on R + x Uq satisfying go(0, u) = and 
(7o(^ ; + x > for x > and u G [To; 

• (x, u) I—)- <7i(x, n) is a Borel function on R + x U\ satisfying g\(x, u) + x > 
for a; > and u G £7i. 

Let {VF(ds, du)} be a white noise on (0,oo) x E with intensity ds7r(dz). Let 
{^(ds, du)} and {Ni(ds,du)} be Poisson random measures on (0,oo) x Uq 
and (0,oo) x U\ with intensities ds no(du) and ds/j,i(du), respectively. Sup- 
pose that {W(ds, du)}, {No(ds,du)} and {Ni(ds,du)} are defined on some 
complete probability space (fl, JP, P) and are independent of each other. Let 
{No(ds, du)} denote the compensated measure of {No(ds, du)}. A nonnega- 
tive cadlag process {x(t) :t > 0} is called a solution of 



if it satisfies the stochastic equation almost surely for every t > 0. We 
say {x(t):t > 0} is a strong solution if, in addition, it is adapted to the 
augmented natural filtration generated by {W(ds,du)}, {No(ds,du)} and 
{Ni(ds, du)} [see, e.g., Situ (2005), page 76]. Since x(s-) ^ x(s) for at most 
countably many s > 0, we can also use x(s) instead of x(s— ) in the integrals 
with respect to W(ds,du) and ds on the right-hand side of (2.1). For the 
convenience of the statements of the results, we write b(x) = bi(x) — b2(x), 
where x \— > b\(x) is continuous, and x \— > 62 (x) is continuous and nondecreas- 
ing. Let us formulate the following conditions: 

(2. a) there is a constant K > so that 



for every x > 0; 

(2.b) there is a nondecreasing function x 1— > L(x) on ]R + and a Borel 




(2.1) 







for every x > 0; 



STOCHASTIC EQUATIONS, FLOWS AND PROCESSES 7 

(2.c) for each m > 1 there is a nondecreasing concave function z \- > r m {z) 
on M_|_ such that J , r m (z)~ 1 dz = oo and 



- + / \g 1 (x,u)-g 1 (y,u)\p,i(du)<r m (\x-y\) 

JUi 

for every < x, y < m; 

(2.d) for each m > 1 there is a nonnegative nondecreasing function z i— > 
Pm(z) on R + so that f Q , p m {z)~ 2 dz = oo, 



|o-(a;,u) - cr(y, u)\ 2 Tr{du) < p m (\x - y\) 2 



and 



/ 



p (du) f 1 l ^ x ^^ 2{l ~^ l {M^v,u)\<n} dt < , , 
' Jo Pm(\(x-y) + tl (x,y,u)\) 2 

for every n > 1 and < x,y < m, where £o0e> y, u) = go(x, u) — go(y, u) and 
c(m,n) > is a constant. 

Theorem 2.1. Suppose that (o~,b,go,gi) are admissible parameters sat- 
isfying conditions (2.a)-(2.d). Then the pathwise uniqueness of solutions 
holds for (2.1). 

Proof. We first fix the integer m > 1. Let ao = 1 and choose at — > 
decreasingly so that J a k ~ 1 Pm{z)~ 2 dz = k for k > 1. Let x i— )■ ipk{x) be a non- 
negative continuous function on R which has support in (afc,afc-i) and satis- 
fies J^' 1 tpk(x) dx = 1 and < ipk{x) < 2k~ l p m {x)~ 2 for a*. < x < afc_i. For 
each k > 1 we define the nonnegative and twice continuously differentiable 
function 

(2.2) 4>k{z)= [ ' dy [ V Mx)dx, z£R. 



It is easy to see that 4>k(z) — > \z\ nondecreasingly as k — > oo and < 4>' k {z) < 1 
for z > and — 1 < <j)'u(z) < for z < 0. By condition (2.d) and the choice of 

Xh-^ fc (x), 

< ^fc(k - 2/|)Pm(k - y|) 2 < - 

for < x,y < m. Then the left-hand side tends to zero uniformly in < 
x, y < m as k — > oo. For /i,(sR, by Taylor's expansion we have 

D h M0= f h 2 <t>'l(( + th)(i-t)dt= f h 2 M\C + th\)(i-t)dt. 

Jo Jo 
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It follows that 

(2.4) D h MO <\f Q h 2 Pm (\( + th\y 2 (l - t) dt. 
Observe also that 

(2.5) D h MO = Afcfo(C) - 4'kiOh < 2\h\. 

For < x, y < m and n > 1 we can use (2.4) and (2.5) to get 

D l {x,y,u)<f>k{x ~ y)va{du) 

Uo 



(2.6) 



kJuo Jo Pm{\{x-y) + tl (x,y,u)\) 



+ 2 1 \l (x,y,u)\lj\ lo r X) y iU )\ >n \iJ,o(du) 



<jc(m,n) +4 / go(m,u)lfo( m ,u)>n/2}t4)(du) 
k Ju 



By conditions (2.b), (2.d) one sees the right-hand side tends to zero uni- 
formly in < x,y < m as k — > oo. Then the pathwise uniqueness for (2.1) 
follows by a trivial modification of Theorem 3.1 in Fu and Li (2010). □ 

The key difference between the above theorem and Theorems 3.2 and 3.3 
of Fu and Li (2010) is that here we do not assume x i— > go(%, u) is nondecreas- 
ing. This is essential for the applications to stochastic equations like (1.6). 

Theorem 2.2. Let (a, b',gQ,g[) and (cr,b" ,go,g'{) be two sets of admis- 
sible parameters satisfying conditions (2.a)-(2.d). In addition, assume that: 

(i) for every u E U\ , x i— )• x + g[(x , u) or x i— )• x + g'{(x,u) is nondecreas- 

ing; 

(ii) b'{x) < b"(x) and g[(x,u) < g'{(x,u) for every x > and u G U\. 

Suppose that {x'(t) :t>0} is a solution of (2.1) with (b,gi) = (b',g[), and 
{x"(t) :t > 0} is a solution of the equation with (b,g\) = (b",g'{). If x'(0) < 
x"(0), then P{x'(t) < x"(t) for allt>0} = l. 

Proof. Let ((t) = x'(t) - x"(t) for t > 0. Let x H> ^(x) be defined as 
in the proof of Theorem 2.1. Instead of (2.2), for each k > 1 we now define 

(2.7) Mz)= [ dy[ V Mx)dx, z£R. 



JO JO 

Then cj>k{z) — > z + : = V z nondecreasingly as k — > oo. Let 

l (t,u) =g {x'{t),u) - g (x"(t),u), t>0,ueU , 
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and 

h(t, u) = g[ (x'{t),u) - g1(x"{t),u), t>0,u£U\. 

For ((s— ) < we have 4>k(C( s ~)) = ^fc(C( s— )) = 0- Since x >->■ x + f(x,u) 
is nondecreasing for f = g[ or g", for ((s— ) = x'(s-) — x"(s—) < we also 
have 

C(s-) + h(s-,u) = x'(s-) - x"(s-) + g ' 1 ( x '( s -),u) - g'l(x"(s-),u) 

< x'(s-) - x"(s-) + f{x'{s-),u) - f(x"{s-),u) < 0. 

The latter implies 

A Jl(s _ tt) 0fc(C(a-)) = MC(s-) + h(s-,u)) - MC(s-)) = 0. 
By Ito's formula we have 

MC{t)) = MC(0)) + \f dsJjUC(s-))[a(x'(s-),u) 

-a{x"{s-),u)] 2 TT(du) 

/V*(C(*-))P<V('-)) 
Jo 

-b"(x"(s-))]l {as _ )>0} ds 

+ / ds A Zl(s _ iU)( /) fc (C(s-))l {c{s _ )>0 }/_fi(du) 
JO JUi 

+ [ ds [ D lo{s _, u)( j) k {({s-))Mdu) + M m {t), 

Jo Ju 

where 

M m (t)= f [ cj>' k (as-))ia(x'(s-),u) 

Jo J E 

-a{x"{s-),u)]W(ds,du) 

+ [ [ A ll{s _ )U)( /) k {C{s-))Ni(ds,du) 
Jo JUi 

+ f [ A lo(s _ u) (j) k (C(s-))No(ds,du). 

Jo Ju 

Let T m = inf{t > : x'{t) > m or x"(t) > m} for m > 1. Under conditions (2.b), 
(2.c) it is easy to show that {M m (t Ar m )} is a martingale. Recall that b'(x) < 
b"(x) and b'{x) = b'^x) — b' 2 (x) for a nondecreasing function x i— > 62 Then 



+ 



(2.8) 
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under the restriction C( s— ) > we have 



<f>' k (((s-))[b'(x'(s-))-b"(x"(s-))} 

< (f> ' k (C(s-))[b'(x'(s-))-b'(x"(s-))} 

< ^ids-mix' (s-^-b^is-))] 

<\b' 1 (x'(s-))-b' 1 (x"(s-))\ 



and 



^i( 8 -,tt)^fc(C(«-)) 

= MC(s-) + g[(x\s-),u)-g'{(x"( S -), U )) 
< MC(s-) + g[(x'( S -),u)-g[(x"( S -),u)) 
<\g' 1 ( x '(s-),n)-g[(x"(s-),n)\. 



The estimates (2.3) and (2.6) are still valid. If x'(0) < x"(0), we can take 
the expectation in (2.8) and let k — > oo to get 

r rt/\T m "1 



where the second inequality holds by the concaveness of z i— )■ r m (z). Then 
E[C(t A r m ) + ] = for all t > 0. Since r m — >• oo as m — >• oo, we get the desired 
comparison property. □ 

We say the comparison property of solutions holds for (2.1) if for any two 
solutions {xi(t) :t > 0} and {^(i) 0} satisfying xi(0) < X2(0) we have 
P{xi(t) < X2(t) for all t > 0} = 1. From Theorem 2.2 we get the following: 

Theorem 2.3. Let (o~,b,go,gi) be admissible parameters satisfying con- 
ditions (2.a)-(2.d). In addition, assume that for every uEUi the function 
x i— > x + gi(x, u) is nondecreasing. Then the comparison property holds for 
the solutions of (2.1). 

The monotonicity assumption on the function x \— > x + gi(x, u) in The- 
orem 2.3 is natural. To see this, suppose that {x\(t)} and {^(t)} are two 
solutions of (2.1) and {(si,Ui):i > 1} is the set of atoms of {Ni(ds,du)}. 
The assumption guarantees that Xi(sj— ) < X2(sj— ) implies 






xi(si) = xi(si-) + gi(x 1 (s i -),u i ) 

< x 2 (si-) + gi(x 2 (si-),Ui) =x 2 (si). 
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A similar explanation can be given to Theorem 2.2. In some applications 
the kernel x i— >■ go(x,u) may be nondecreasing. When this is true, we can 
replace (2.d) by the following simpler condition: 

(2.e) For each u€Uq the function x \— > go(x, u) is nondecreasing, and for 
each m > 1 there is a nonnegative and nondecreasing function z i— >■ p m {z) 
on R + so that J 0+ p m (z)~ 2 dz = oo and 



cr(cc,u) -cr(y,u)| 7r(dii) + / |Z (^, 2/, ^)| A |Z (^,y,«)| Mo(d«) 
<Pm.{\x-y\) 2 

for all < x, y < m, where lo(x, y, u) = go(x, u) — go{y, u). 

Proposition 2.4. Let (a,b,go,gi) be admissible parameters. If (2.e) 
holds, then (2.d) holds. 

Proof. Since x i-> go(x,u) is nondecreasing, it is not hard to see \(x — 
y) + Uq(x, y, u)\ >\x — y\. By condition (2.e) and the monotonicity of z h- > 
p(z) we have 

Cdtt 

7o M ^m(K^-y) + «o(a;,y,«)|) 2 

. A 1 f [Mx,y,u)\Al (x,y,u) 2 } 

<n dt po(du)<n. 

Jo Juo Pm{\x-y\Y 

Then condition (2.d) is satisfied. □ 

Theorem 2.5. Suppose that (o~,b,go,gi) are admissible parameters sat- 
isfying conditions (2. a), (2.c), (2.e). Then there is a unique strong solution 
to (2.1). 

Proof. We first note that (2.b) follows from (2.e). By Proposition 2.4, 
we also have (2.d) from (2.e). Let {V n } be a nondecreasing sequence of 
Borel subsets of Uo so that lj^=i = Uo and Ho(V n ) < oo for every n > 1. 
For m,n>l one can use (2.e) to see 

x^p m (x) := / [g (x,u) - g (x,u) Am]po(du) 
JUo 

and 

x H> j m ,n( x ) '■= / [9o(x, u) A m]no{du) 
Jv n 

are continuous nondecreasing functions. By the results for continuous-type 
stochastic equations as in Ikeda and Watanabe [(1989), page 169] one can 
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show there is a nonnegative weak solution to 

x(t) = x(0)+ [ [ a(x(s) Am,u)W(ds,du) 

JO J E 

+ / b m (x(s) A m) ds - 7 m ,„(x(s) Am)ds, 
Jo Jo 

where b m (x) = b(x) — (3 m (x). The pathwise uniqueness holds for the above 
equation by Theorem 2.1. Then it has a unique strong solution. Let {VF n } be 
a nondecreasing sequence of Borel subsets of U\ so that IJ^Li W n = U\ and 
A*i (Wn) < 00 f° r every n > 1. Following the proof of Proposition 2.2 of Fu and 
Li (2010) one can show there is a unique strong solution {x m n (t) : t > 0} to 

x(t)=x(0)+ I I a(x(s-) Am,u)W(ds,du) 
Jo J E 

+ I b m (x(s-) Am)ds - I "f m ,n(x(s) A m) ds 
Jo Jo 

/ [9o(x(s— ) A m, u) A m]No(ds, du) 
'0 Jv n 

+ / / [di(x(s— ) A m, u) A m]Ni(ds, du). 
Jo Jw n 

We can rewrite the above equation into 

x(t)=x(0)+ [ [ a(x(s-) Am,u)W(ds,du) 
Jo J E 

+ / b m (x{s—) Am) ds 
Jo 

+ / / [go(x(s— ) A m, u) A m]No(ds, du) 

Jo Jv n 

+ / / [si(x(s—) A m, u) A m]Ni (ds, du). 
Jo Jw„ 

As in the proof of Lemma 4.3 of Fu and Li (2010) one can see the sequence 
{x m ,n(t) - t > 0}, n = 1,2, . .. , is tight in D([0, oo),R+). Following the proof 
of Theorem 4.4 of Fu and Li (2010) it is easy to show that any weak limit 
point {x m (t) : t > 0} of the sequence is a nonnegative weak solution to 

x(t) = x(0)+ / / a(x(s-) Am,u)W(ds,du) 
Jo J E 

+ / b m (x(s-) Am)ds 
Jo 
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(2.9) 

+ / / [do( x {s— ) A to, u) A m]No(ds, du) 
Jo Ju 

+ / / [gi(x(s—) Am,u) Am]Ni(ds,du). 
Jo JlJx 

By Theorem 2.1 the pathwise uniqueness holds for (2.9), so the equation has 
a unique strong solution [see, e.g., Situ (2005), page 104]. Then the result 
follows by a simple modification of the proof of Proposition 2.4 of Fu and Li 
(2010). See El Karoui and Meleard (1990) and Kurtz (2007, 2010) for the 
general theory of stochastic equations driven by white noises and Poisson 
random measures. □ 

3. Stochastic flows of CBI-processes. In this section, we give the con- 
structions and characterizations of the flow of CBI-processes and the asso- 
ciated immigration superprocess. Suppose that a > and b are constants, 
and (u A u 2 )m(du) is a finite measure on (0, oo). Let (j> be a function given 
by 

l f°° 

(3.1) (j)(z) = bz + -a 2 z 2 + / (e~ zu - 1 + zu)m(du), z>0. 

A Markov process with state space ffi+ := [0, oo) is called a CB-process with 
branching mechanism eft if it has transition semigroup (pt)t>o given by 

(3.2) f e- x »p t (x,dy) = e- XVt W, A > 0, 

where (t, A) i— > ft (A) is the unique locally bounded nonnegative solution of 

j t v t (X) = -<f>(vtW), v (X) = X, t>0. 
Given any (3 > we can also define a transition semigroup (qt)t>o on R + by 

(3.3) / e- Xv q t (x,dy) = exd-xv t (\)- f Pv s (X)ds 



A nonnegative real- valued Markov process with transition semigroup (qt)t>o 
is called a CBI-process with branching mechanism (j) and immigration rate (5. 
It is easy to see that both (jpt)t>o an d {qt)t>o are Feller semigroups. See, for 
example, Kawazu and Watanabe (1971) and Li (2011), Chapter 3. 

Let {W(ds, du)} be a white noise on (0, oo) 2 based on the Lebesgue mea- 
sure dsdu, and let { N(ds, dz, du)} be Poisson random measure on (0,oo) 3 
with intensity dsm(dz) du. Let {N(ds,dz,du)} be the compensated measure 
of {N(ds, dz, du)}. 
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Theorem 3.1. There is a unique nonnegative strong solution of the 
stochastic equation 

ft r-Y a - ft 

Y t = Y + a / W(ds, du) + / {fi - bY s _) ds 
Jo Jo Jo 

ft foo fY s - 

+ / zN(ds,dz,du). 

Jo Jo Jo 

Moreover, the solution {Yt : t > 0} is a CBI-process with branching mecha- 
nism <j> and immigration rate (3. 

Proof. The existence and uniqueness of the strong solution follows by 
an application of Theorem 2.5 [see also Dawson and Li (2006)]. Using Ito's 
formula one can see that {Yt(v) :t > 0} solves the martingale problem asso- 
ciated with the generator L defined by 

1 

(3.4) Lf(x) = -a 2 xf"{x) + (f3-bx)f'(x)+x DJ(x)m(dz). 



Then it is a CBI-process with branching mechanism <fi and immigration 
rate j3 [see Kawazu and Watanabe (1971) and Li (2011), Section 9.5]. □ 

Let v I—)- j(v) be a nonnegative and nondecreasing continuous function on 
[0,oo). We denote by j(dv) the Radon measure on [0, oo) so that 7([0,w]) = 
"f(v) for v > 0. By Theorem 3.1 for each v > there is a pathwise unique 
nonnegative solution Y(v) = {Yt(v) :t > 0} to the stochastic equation 

ft fY s -(v) ft 

Y t (v)=v + a / W(ds,du)+ [j(v)-bY s „(v)]ds 
(3.5) J ° J ° J ° 

ft foo fY s -(v) 

+ 11 \ zN(ds,dz,du). 



JO Jo 



Theorem 3.2. For any V2>vi>0 we have P{Y t (v2) > Y t (v\) for all 
t > 0} = 1 and {Y t {v2) — Yt(v\) :t > 0} is a CBI-process with branching mech- 
anism <j) and immigration rate (3 := 7(1*2 ) — 7(^1) ^ 0. 

Proof. The comparison property follows by applying Theorem 2.2 and 
Proposition 2.4 to (3.5). Let Z t = Y t {v2) — Y t {v\) for t > 0. From (3.5) we have 

ft fY s -(v2) ft 

Z t = v 2 -v 1 + a / W(ds,du)+ / (p-bZ s -)ds 

Jo Jo 

ft foo fY s -(v2) 

+ / / / zN(ds,dz,du) 
Jo Jo Jy s -{vi) 



(3.6) 
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rt rZ 3 - rt 

= v 2 -vi+a / Wi(ds,du)+ / (p-bZ s J)ds 
Jo Jo Jo 

rt roo rZ s - 

+ / zNi(ds,dz,du)i 

Jo Jo Jo 



where 



Wi {ds, du) = W(ds, y s _ («i ) + du) 

is a white noise with intensity ds du, and 

Ni(ds, dz, du) = N(ds, dz, Y s -(vi) + du) 

is a Poisson random measure with intensity ds m{dz) du. That shows {Zt : t > 
0} is a weak solution of (3.5). Then it a CBI-process with branching mech- 
anism <j) and immigration rate (3. □ 

Theorem 3.3. Let v 2 >v 1 >u 2 >u 1 > 0. Then {Y t (u 2 ) -Y t (u x ):t> 0} 
and {Yt(v 2 ) — Yt(v\) : t > 0} are independent CBI-processes with immigration 
rates a := 7(1*2) — 7(^1) and f3 := 7(^2) — 7(1*1), respectively. 

Proof. Let L a and Lp denote the generators of the CBI-processes with 
immigration rates a and /?, respectively. Let X t = Y t (u 2 ) — Y t {u\) and Z t = 
^(^2) — Yt{v\). For any G £ C 2 (IR^_) one can use Ito's formula to show 

rt 



(3.7) 



G{X t ,Z t ) = G(X ,Z )+ [ L a G(X s ,Z s )ds 

Jo 

+ / LpG(X s , Z s ) ds + local mart., 
./o 



where L a and Lp act on the first and second coordinates of G, respectively. 
Then {Xt : t > 0} and {Z( : t > 0} are independent CBI-processes with immi- 
gration rates a and /3, respectively. □ 

Proposition 3.4. There is a locally bounded nonnegative function 1 1— > 
C(t) on [0,oo) so that 

E{ sup [Y s (v 2 )-Y s (v 1 )]\<C(t){(v 2 -vi) + ty(v 2 )-~,(v 1 )] 

l 0<s<t > 

(3.8) 



+ yjv 2 - Vi + yj 7(^2) - 7(^l)} 

for t > and v 2 >v\> 0. 

Proof. Let Z t = Y t (v 2 ) — Y t (v\) for t > 0. Taking the expectation in (3.6) 
we have 

E(Z t ) = («2 - vx) + t[j(v 2 ) - 7 ( V1 )] - 6 / E(Z S ) ds. 

j 
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Solving the above integral equation gives 

(3.9) E(Z t ) = (v 2 - Vl )e- bt + [7(^2) - 7(«i)]& _1 (l " 

with 6 _1 (1 — e~ bt ) = t for b = by convention. By (3.6) and Doob's martin- 
gale inequality, 



E<| sup Z s }< («2 - vi) + 2aE 1/2 



t pY a -{v 2 ) \ 2 

/ W(ds,du) 

JYs-iin) 



+ 



+ 



2E 1 / 2 



{[7(«2)-7(«i)] + |6|E(^.)}ds 
t /.1 /•n_(« 2 ) 



+ E 



JO JY B -(vi) 
t poo pY s -(v 2 ) 

3 Jx 7y s _(«i) 



zN{ds, dz, du) 
zN{ds, dz, du) 



< (v2 ~ vi) + t[j(v 2 ) - 7(^1)] + 2ct 

1/2 



+ 2 



z v{dz) 



E(Z s )ds 



E(Z S ) ds 

1/2 



1/2 



+ 



zv(dz) 



E(Z s )ds. 



Then (3.8) follows by (3.9). □ 



Suppose that (E, p) is a complete metric space. Let F be a subset of [0, 00) 
such that £ F and let 1 1— )• be a path from i 7 to E. For any e > the 
number of e-oscillations of this path on F is defined as 

pi(e) := sup{n > : there are = to < h < • ■ • < t n G F 

so that p(x(ti-i),x(ti)) > e for all 1 < i < n}. 

If F is dense in [0, 00), it is simple to show the limits y(t) := ]im.FB8-rt+ x(s) 
exist for all t > and constitute a cadlag path t >— > y(t) on [0, 00) if and only 
if t 1 — y x(t) has at most a finite number of e-oscillations on Ff] [0, T] for every 
e > and T > 0. 



Lemma 3.5. Suppose that (J),£f,5ft,P) is a filtered probability space and 
{Xt :t > 0} is o (^t)-Markov process with state space {E,S) and transition 
semigroup {P s ,t)t>s- Suppose that p is a complete metric on E so that: 

(i) for £ > and < s,t < u we have {u £ f2 : p(X s (uj), X t (uj)) < e} S 
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(ii) for e > and x G E we have U £ (x) := {y G E : p(x, y) < e} G £ and 
(3.10) a £ (h):= sup sup P s , t (x, U £ {x) c ) -> (/i->0). 

0<t-s<hx€E 

Then {Xt : t > 0} /ias a p-cadlag modification. 

Proof. Let F = {0,ri,r2, • . .} be a countable dense subset of [0,oo) 
and let F n = {0, r±, . . . , r n }. For e > and a > let z^ a (e) and z-^(e) denote, 
respectively, the numbers of e-oscillations of t \— > Xt on F n [0, a] and F n n 
[0, a]. Then v^(e) — > v a {e) increasingly as n — > oo. Let r^(0) = and for k > 
define 

+ 1) = mm{t G F n n «(/c), oo) : p(X Tk{k) , X t ) > e}, 

if r^(A;) < oo and r^(A; + 1) = oo if r^{k) = oo. Since F n is discrete, for any 
a > we have 

+ 1) < a} = |J = S } n {p(X s ,X t ) > e }). 

s<tGF n n[0,a] 

Using property (i) and the above relation it is easy to see successively that 
each T^(k) is a stopping time. As in the proof of Lemma 9.1 of Wentzell 
[(1981), page 168] one can prove P{r^(l) < h} < 2a £ / 2 {h) for e > and 
h > 0. Since the strong Markov property of {Xt : i > 0} holds at the discrete 
stopping times r^(fc), A; = 1, 2, ... , one can inductively show 

P{^(2e) > A;} < P«(fe) < /i} < P^fa)]*- 

It follows that 

P{v h (2s) >k}= lim P{^(2e) > A;} < [2a e/2 (/0]* 

Choosing sufficiently small h = h(e) G Fn (0, oo) so that a £ /2(h) < 1/2 and 
letting k — > oo we get P{^ ,ft, (2e) < oo} = 1. By repeating the above pro- 
cedure successively on the intervals [h,2h], [2h, 3/i], . . . we get P{z^ a (2e) < 
oo} = 1 for every a > 0. Let fti = fl^ii^lVm) < oo}. Then ^ G Sf and 
P(Oi) = 1. Moreover, for £j G we can define a p-cadlag path t *— > Yt(u}) 
on [0,oo) by Y t (oj) := limF3 S -*t+X s (oj). Take xq G F and define lt(w) = xo 
for t > and w G \ Sli- By (3.10) one can see t >— > Xt is right continuous 
in probability, so Yt = Xt a.s. for every t > 0. Then {Y^ : t > 0} is a p-cadlag 
modification of : t > 0}. □ 

Let D[0,oo) be the space of nonnegative cadlag functions on [0,oo), and 
let 3$(D[Q,oo)) be its Borel a-algebra generated by the Skorokhod topology. 
Theorems 3.2 and 3.3 imply that {Y(v) : v > 0} is a nondecreasing process 
in (D[0, oo), &(D[0, oo))) with independent increments. Let p be the metric 
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on D[0,oo) denned by 

/•oo 

(3.11) p(C, C) = / e~ l sup - CO0I A 1) (it. 

JO 0<s<t 

This metric corresponds to the topology of local uniform convergence, which 
is strictly stronger than the Skorokhod topology. 

Theorem 3.6. The path-valued process {Y(v) : v > 0} has a p-cadlag mo- 
dification. Consequently, there is a version of the solution flow {Yt{v) : t > 0, 
v > 0} o/ (5.5) witt the following properties: 

(i) for each v > 0, f i- >• lt(u) is a cadlag process on [0,oo) and solves (3.5); 

(ii) for each t>0,v*-^Yt(v) is a nonnegative and nondecreasing cadlag 
process on [0, oo). 

Proof. Step 1. For any T > let D[0,T] be the space of nonnegative 
cadlag functions on [0, T], and let &(D[0,T]) be its a-algebra generated 
by the Skorokhod topology. For v > let Y T (v) = {Y t (v) : < t < T}. The- 
orem 3.3 implies that {Y T (v):v > 0} is a process in (D[0,T],&(D[0,T])) 
with independent increments. 

Step 2. Let = {T,ri,r<2, . . .} be a countable dense subset of [0,T]. We 
consider the metric px on Z)[0,T] defined by 

M£,0= sup - C(s)| = sup - C(«)l- 

0<s<T r£F T 

For any e > and £ G -D[0, T] we have 

^ e (e):={C€D[0,r|:/> r (^ > C)<e} 

= f| {CeZ>[0,r]:|er-Cr|<e}. 

r£F T 

Then the above set belongs to &(D[Q, T]) [see, e.g., Ethier and Kurtz (1986), 
page 127]. It follows that 

oo 

U E (0 := {C € D[0, T] : p T (H, () < e} = (J ^ e _ 1/n (0 

71=1 

also belongs to ^(Z?[0,T]). 

Step 3. Let (^J) v >o be the natural filtration of {y T (t>) : v > 0}. For any 
e > and < s, i < v we have 

,»r(Y T (s),y T (t)) = sup \Y r (s)-Y r (t)\. 

r<=F T 

Then one can show {w G : p T (Y T (uj, s), Y t (lo, t)) < e} G 
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Step 4. Let (P^ v ) v >u denote the transition semigroup of {Y T (v) :v > 0}. 
By Proposition 3.4 for e > and £ £ -D[0, oo) we have 



PuA^U £ (O c )=~P{ sup [Y s (v)-Y s (u)]>e\ 
Ks^eI sup [Y s (v)-Y s (u)]\ 

l 0<s<T ' 

<e- 1 C(t){(v-u) + h(v)->y(u)] 



+ y/v — u + yj(v) - j(u)}. 

Since v h- > 7(11) is uniformly continuous on each bounded interval, Lemma 3.5 
implies that {Y T (v) :v > 0} has a p^-cadlag modification. That implies the 
existence of a p-cadlag modification of {Y(v) :v > 0}. □ 

In the situation of Theorem 3.6 we call the solution {Yt(v) :t > 0,v > 0} 
of (3.5) a flow of CBI-processes. Let F[0, 00) be the set of nonnegative and 
nondecreasing cadlag functions on [0, 00). Given a finite stopping time r and 
a function fi £ F[0, 00) let {Y^ t (v) :t > 0} be the solution of 

Y£ t (v) = u(v) + J J W(ds, du) 

(3.12) + r +t Mv)-bY»(v)}ds 



pr+t roo rYf s _(v) 

+ / / zN(ds,dz,du) 

Jt Jo Jo 

and write simply {Y^(v):t > 0} instead of {YQ t (v):t > 0}. The pathwise 
uniqueness for the above equation follows from that of (3.5) since {VF(r + 
ds, du)} is a white noise based on ds dz, and {N(t + ds,dz, du)} is a Poisson 
random measure with intensity ds m(dz) du. Let G T j be the random operator 
on .F[0,oo) that maps u to Y^ t . 

Theorem 3.7. For any finite stopping time r we have P{Y^ +t = G Ti tYf 
for all t >0} = 1. 

Proof. By the sample path regularity of (t, v) \— > Yt(v) we only need 
to show P{Yf +t (v) = G Tj tYr(v)} = 1 for every t > and v > 0. By (3.5) we 
have 

Yf +t (v) = Y?{v) + a / W(ds, du) 



r +t [j(v)-bY^(v)]ds 
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i-r+t r-co rY^_{v) 

+ / / / zN(ds,dz,du). 
Jt Jo Jo 

By the pathwise uniqueness for (3.12) we get the desired result. □ 

For any Radon measure n(dv) on [0, oo) with distribution function v\-t fJ-(v), 
the random function v t— > Y^{y) induces a random Radon measure Y^dv) 
on [0, oo) so that Y t ^([0, v]) = Y^iy) for v > 0. We shall give some character- 
izations of the measure- valued process {Y/ 1 : t > 0}. 

For simplicity, we fix a constant a > and consider the restrictions of fJ,(dv), 
"f(dv) and {Y t ^:t > 0} to [0, a] without changing the notation. Let us con- 
sider the step function 

n 

(3.13) f(x) =c l{ }{x) + J2 c i l (a % - uai ](.x), xG[0,a], 

i=l 

where {cq, c%, . . . , c n } C R and {0 = ao < a\ < • • • < a n = a} is a partition of 
[0,a]. For this function we have 

n 

(3.14) (YfJ) = c o Y t ^0) + ^c i [Yr(a i )-Y t ^ 1 )]. 

i=l 

From (3.12) and (3.14) it is simple to see 

rt roo 

(Yf,f) = ( f iJ) + a / rf_(u)W(da,du) 
Jo Jo 

(3.15) + f[tn,f)-bQT?_,f)]da 

Jo 

rt roo roc 

+ / / / zgg_(u)N(ds,dz,du), 
Jo Jo Jo 

where 

n 

(3.16) g%(u) = C l{ n <yM (0) } + c i 1 W(a i _- l )<u<Yj > (g i )}- 

i=l 

Proposition 3.8. For any t>0 and f £ B[0,a] we have 

(3.17) E[(Yf , /)] = /)e~ w + ( 7 , /)6" 1 (1 - e~ bt ) 
with b~ l (l — e~ ht ) = t for b = by convention. 

Proof. We first consider the step function (3.13). By taking the expec- 
tation in (3.15) we obtain 

E[(Yf, /)] = (jm, f) + i< 7 , f)~bf E[(Y/\ /)] ds. 

Jo 
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The above integral equation has the unique solution given by (3.17). For 
a general function f £ B[0, a] we get (3.17) by a monotone class argument. 

□ 



Theorem 3.9. The measure-valued process {Y^" : t > 0} is a cddldg strong 
Markov process in M[0,a] with Yq = jjL. 

Proof. In view of (3.14), the process 1 1— > (Y^,f) is cadlag for the 
step function (3.13). Since any function in C[0, a] can be approximated by 
a sequence of step functions in the supremum norm, it is easy to conclude 
1 1 — y (Xt i f) i s cadlag for all / G C[0, a]. By Theorem 3.7, for any finite stop- 
ping time r we have Y^ +t = G t jYt almost surely. That clearly implies the 
strong Markov property of {Yj 1 : t > 0}. □ 

Theorem 3.10. For any f G B[0, a] the process {(Yf, f) :t > 0} has 
a cadlag modification. Moreover, there is a locally bounded function 1 1— >■ C(t) 
so that 

(3.18) E[ sup {Y?J)] <C(t)[^,/) + ( 7 ,/) + (/i,/ 2 ) 1/2 + ( 7 ,/ 2 ) 1/2 ] 

l 0<s<t J 

for every i > and f G i?[0,a] + . 

Proof. We first consider a nonnegative step function given by (3.13) 
with constants {cq,ci, . . . ,c n } Cl+. By (3.15) and Doob's martingale in- 
equality, 



E 



sup (Y?,f) 

0<s<t 



<{ f i,f) + 2aB 1 / 2 



t poo 



J 



g^_(u)W(ds,du) 



+ t(7,f) + \b\ I E[(Y»,f)]ds 
Jo 



+ 2E 1/2 
+ E 



t pi roo 



zgg- (u)N(ds, dz, du) 



o Jo Jo 

t POO roo 



Jl JO 



zg^_ (u)N(ds, dz, du) 



^,/)+2aE 1 /2 



ds 



g%(u) 2 du 



+ t(j,f) + \b\ J)]d8 
Jo 
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+ 2E 1/2 
+ E 



ds I z 2 m(dz) I g%(u) 2 du 
o Jo Jo 



t poo 

ds / zm(dz) I g^(u)du 



o J\ 



o 



<(H,f) + 2U V[(Y»J 2 )]ds 



1/2 



r / r 1 


\ V2l 


cr+ / z 2 m(dz) 


L \io 




/>oo 




- / zm(dz) 









In view of (3.17) we get (3.18) for the step function. Now let rj(dv) = /j,(dv) + 
ry(dv) and choose a bounded sequence of step functions {/ n } so that / n — > f 
in L 2 (r^) as n — > oo. By applying (3.18) to the nonnegative step function 
|/n - /ml we get 



E 



sup {Yf, \f n - f m \) < C(t)[( V , \f n -f m \)+ 2(77, |/n - fm?) 1 ' 2 ]. 
0<s<t J 



The right-hand side tends to zero as m,n — > oo. Then there is a cadlag 
process {Yf'(f) : t > 0} so that 



(3.19) 



E 



0, n— >oo. 



sup 

L 0<s<t J 

On the other hand, from (3.17) we have 

E[(lf , \f n -f\)} = (fi, \f n - f\)e~ bt + b-\l - e- bt )( 7 , \f n - /|), 

which tends to zero as n — > oo. Then {Y t ^(f):t > 0} is a modification of 
{(YfJ) -t > 0}. Finally, we get (3.18) for / G B[0,a}+ by using (3.19) and 
the result for step functions. □ 

Theorem 3.11. The process {Yf'-.t > 0} is the unique solution of the 
following martingale problem: for every G G C 2 (M) and f G £>[0, a], 



GW\f)) 



1 



G((n,f)) + ^ I G'\{Y?,f)){Yt,f 2 )ds 



(3.20) 



+ / G'((y/,/))[( 7 ,/)-6(^,/)]d S 
Jo 

+ ds Yj>(dx) / [G{{Y?J)+zf{x)) 

J0 J[0,a] JO 

- G((YfJ)) - zf(x)G'((Y», f))]m(dz) 



+ local mart. 
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Proof. Again we start with the step function (3.13). Using (3.15) and 
Ito's formula, 

i rt roc 

+ [ t G'(<y£_j))[( 7 ,f)-b<y£.j)]d8 

Jo 

ft /*oo /"OO 

+ ds m(dz) / [G((Y»,f) + zg»(u)) 

JO JO J 

-G((y/,/))-G'((y/,/))^(n)](iu 

+ local mart. 
= G((m, /)) + \o 2 f G"((Yf, f))(Yf, f) ds 

+ fG\{Y?J))[( Ci ,f)-b(Y?,f)]ds 
Jo 

+ ds Y?{dx) / [G(<y?,f)+zf(x)) 
Jo Jo Jo 

- G((Yf, /)) - G\{YtJ))zf{x)]m{dz) 

+ local mart. 

That proves (3.20) for step functions. For / € £?[0,a] we get the martingale 
problem using (3.19). The uniqueness of the solution follows from a result 
in Li (2011), Section 9.3. □ 

The solution of the martingale problem (3.20) is the special case of the 
immigration superprocess studied in Li (2011) with trivial spatial motion. 
More precisely, the infinitesimal particles propagate in [0, a] without migra- 
tion. Then for any disjoint bounded Borel subsets B\ and B2 of [0, a], the 
nonnegative real-valued processes {Y/ J, (B\) : t > 0} and \Y^{B2) ■ t > 0} are 
independent. That explains why the restriction of {Y/ 4 : t > 0} to the interval 
[0, a] is still a Markov process. To consider the process of measures on the 
half line [0,oo) we need to introduce a weight function as follows. 

Let h be a strictly positive continuous function on [0, 00) vanishing at 
infinity. Let M/J0, 00) be the space of Radon measures fi on [0,oo) so that 
(n, h) < 00. Let Bh[0, 00) be the set of Borel functions on [0, 00) bounded by 
const • h, and let C/J0, 00) denote its subset of continuous functions. A topo- 
logy on M^[0,oo) can be defined by the convention /i n , — > \i in M/J0, 00) 
if and only if (/i„, /) — > (fi, /) for every / £ C/J0, 00). Suppose that \x G 
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M/j[0, oo) and 7 S M/JO, 00). It is easy to show that {Yf : t > 0} is a cadlag 
strong Markov process in M/JO, 00), and the results of Theorem 3.10 and 
Theorem 3.11 are also true for 5^(0, 00). 

4. Generalized Fleming Viot flows. In this section we give a construc- 
tion of the generalized Fleming-Viot flow as the strong solution of a stochas- 
tic integral equation. Let a > 0, b > and < (3 < 1 be constants, and 
let z 2 u{dz) be a finite measure on (0, 1]. Suppose that {B(ds, du)} is a white 
noise on (0,oo) 2 with intensity dsdu, and {M(ds,dz,du)} is a Poisson ran- 
dom measure on (0,oo) x (0,1] x (0,oo) with intensity dsv(dz) du. Let 

q(x,u) = l {u < 1Ax} - (1 Ax), x > 0,u E (0,1]. 

We first consider the stochastic integral equation 

X t = X + [ [ aq(X s _,u)B(ds,du) + [ b{f3-X s _)ds 
(41) JoJo^ Jo 

+ zq(X s -,u)M(ds,dz,du), 
Jo Jo Jo 

where M(ds,dz,du) denotes the compensated measure of M(ds,dz,du). In 
fact, the compensation in (4.1) can be disregarded as 

f q(X s _,u)du= [ [l {u < Xs _ A1} - (AV Al)]d« = 0. 
Jo Jo 



Theorem 4.1. There is a unique nonnegative strong solution to (4-1)- 

Proof. We first show the pathwise uniqueness for (4.1). Set l(x,y,u) = 
q{x, u) — q(y, u). For x, y > and < z, t < 1 we have 

(x - y) + ztl(x,y,u) 

= [{x - 1 A x) - (y - 1 A y)} + (1 - zt)(l A x - 1 A y) 

+ z t(~L{u<xAl} ~ l{w<2/Al})- 

It is then easy to see 

\(x - y) + ztl(x,y,u)\ > (1 — zt)\l A x — 1 A y\. 
Moreover, we have 

/ l(x, y, u) 2 du = (1 A x - 1 A y) - (1 A x - 1 A y) 2 
Jo 

< |1 Ax — 1 Ay\. 
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Using the above two inequalities, 



(l-t)dt / v{dz) 



z 2 l(x,y,uf 



< 



< 



l v(dz) 



l u(dz) 



\(x 

- 1 



zt 

- 1 



dt 



y) + ztl(x,y,u)\ 
1 l(x,y,u) 2 



du 



o 



l-zt 



■dt < 



I Ax 
l 

2 



lAy| 
z z v{dz). 



■ du 



o 



Then condition (2.d) is satisfied with p{z) = \fz. Other conditions of The- 
orem 2.1 can be checked easily. Then we have the pathwise uniqueness 
for (4.1). To show the existence of the solution, we may assume Xq = v > is 
a deterministic constant. By Theorem 2.5 there a unique nonnegative strong 
solution of (4.1) if the Poisson integral term is removed. Then for each k > 1 
there is a unique nonnegative strong solution to 



Z + 



JO 



aq(Z s ^,u)B(ds, du) 



(4.2) 



+ / b(P-Z s _)ds 



Jl/kJO 



zq(Z s _,u)M(ds, dz, du), 



because the last term on the right-hand side gives at most a finite number 
of jumps on each bounded time interval. Let {Zk(t) :t > 0} be the solution 
of (4.2) with Z k (0) = v. Let T\ = inf{t > : Z k {t) < 1}. On the time interval 
[0,Ti], the stochastic integral terms in (4.2) vanish. Then 1 1— )■ Z k {t) is non- 
increasing on [0, Ti] . By modifying the proof of Proposition 2.1 in Fu and Li 
(2010) one can see Z k (t) < 1 for t > T x . Thus Z k (t) < (Z k (0) VI) = (v V 1) for 
all t > 0. Let {r k } be a bounded sequence of stopping times. Note that the 
last term on the right-hand side of (4.2) can be considered as a stochastic 
integral with respect to the compensated Poisson random measure. Then for 
any t > we have 

B{[Z k ( Tk + t)-Z k (T k )} 2 } 



< 3a 2 B 



3E 



lis / q(Z k (r k + s),u) du 



+ 3b 2 t 2 (v\Jiy 



ds 



z 2 v{dz) \ q(Z k (r k + s),u) du 



< 3t 



a' 



+ tb 2 (v\/l) 2 + 



z 2 v{dz) 
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The right-hand side tends to zero as t — > 0. By a criterion of Aldous (1978), 
the sequence {Zk(t):t > 0} is tight in D([0, oo), M+) [see also Ethier and 
Kurtz (1986), pages 137 and 138]. By a modification of the proof of Theo- 
rem 4.4 in Fu and Li (2010) one sees that any limit point of this sequence 
is a weak solution of (4.1). □ 

Now let v i— > 7(f) be a nondecreasing continuous function on [0, 1] so 
that < 7(1;) < 1 for all < v < 1. We denote by j(dv ) the sub-probability 
measure on [0, 1] so that 7([0,v]) = 7(f) for < v < 1. By Theorem 4.1 for 
each v > there is a pathwise unique nonnegative solution {X t (v ) : t > 0} to 
the equation 



It is not hard to see that < v < 1 implies P{0 < X t (v) < 1 for all t > 0} = 1. 
The compensation for the Poisson random measure can be disregarded, so 
this equation just coincides with (1.6). By Theorem 2.2 for any < v\ < 
«2 < 1 we have 

P-PQ(ui) < X t (v 2 ) for all t > 0} = 1. 
Therefore {^(u) : < w < 1} is a nondecreasing path- valued process in L>[0, 00). 

Proposition 4.2. There is a locally bounded nonnegative function t*-> 
C(t) on [0,oo) so that 



+ y/v 2 -vi + y/l{v 2 ) - 7(^1)} 

for t>0 and < v\ < v 2 < 1 . 

Proof. Let Z t = X t (v 2 ) - X t (v\) for t > 0. From (4.3) we have 



(4.3) 





(4.4) 



<C7(t){(«2-^) + [7(«2)-7(ui)] 



(4.5) 
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where Y s (u) = L{x s (v 1 )<u<x s (v 2 )}- Taking the expectation in (4.5) and solving 
a deterministic integral equation one can show 

(4.6) E[Z t ] = (v 2 - Vl )e- bt + [ 7 fa) - 7(«l)](l - e~ bt ). 

By (4.5) and Doob's martingale inequality, 

E^sup^ Z s } < {v 2 - v x ) + 2CTE 1 / 2 1 Qf jT 1 [F s _ (n) - Z s _].B(d«, | 

+ /"*6{[7(w2)-7(«i)] + E[Z a ]}da 
■/ o 

+ 2E 1 / 2 jYjf jf J z[Y s .(u) - Z s -)M(ds,dz,du)\ \ 
= (v 2 - v 1 ) + 2aE 1 / 2 ^ds J\y s (u) - Z s 

JO 



' du 



z 2 v{dz) I [Y s (u) -Z s ] 2 du 



where 



L 



i 

\Y s {u) - Z s f du = Z s (l - Z a ) < Z s . 



'o 

Then we have (4.4) by (4.6). □ 

Recall that D[0, oo) is the space of nonnegative cadlag functions on [0, oo) 
endowed with the Borel <r-algebra generated by the Skorokhod topology. 
Let p be the metric on Z)[0, oo) defined by (3.11). 

Theorem 4.3. The path-valued process {X{v) :0 < v < 1} is a Markov 
process in D[0, oo). 

Proof. Let < v < 1, and let r n = mf{t > : -Xt(u) < 1/n} for n > 1. In 
view of (4.3), we have X t {v) = if = 0. Then r n -»■ := inf{t > : 

Xt (v) = 0} as n — > oo. For any p G [0, v) the comparison property and path- 
wise uniqueness for (4.3) imply X t (p) = X t (v) for t > t^. Let Z n {t) = 
XtArM^XtATnip) for t > 0. By (4.3) and Ito's formula, 

p f tATn f 1 a 
Z n (t) = -+ / ———[l {u < Xs{p)} -X s _{p)]B(ds,du) 

v JO JO ^s-\V) 

^ f 1 aX s ^(p) 

'0 
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+ 



+ 



tAT n 



bX a -(v)- 1 [y(p)-7(v)X s -(v)- 1 X s -(p)]d8 



tAT n 



ds 



tATn 



ds 



1 a 2 X s _(p) 2 
1 a 2 



o *^) 2 



[]><x s _(p)} -^s-(p)] 



+ 
P 

+ 
+ 



x [l{u<X s -(v)} ~ X s -{v)]du 

tATn f l f ^X s _(p)(l-z)+Zl {u < Xs _ m X S _( P ) 




+ 



J0 [X s -(v)(l- z) + zl {u < Xs _( v )} X s-( V ). 

x M(ds,dz,du) 

tATn fX s -(v) 

x B(ds, du) 



o J o 



6X s _( V )- 1 [ 7 (p)-7( 7 ;)X s _( ? ;)- 1 X s _(p)]d S 



o 



tATn rl rX s -(v) 

Jo Jo 



Xs-(p)(l -z) + zl {u < Xs _(p)} X s ^(p) 



X,-(v)(l-z) + z 
x M(ds, dz, du), 



X s (v) 



where the two terms involving a 2 counteract each other. Observe also that 
the last integral does not change if we replace M(ds, dz, du) by the compen- 
sated measure M(ds,dz,du). Then we get the equation 



2n(*) = - + 

V Jo 



tAT n fX s -(v) 



aX s _(v) 1 [l{ u <x,_(«)Z tt («-)} - Z n{s-)\ 



x B(ds, du) 



(4.7) 



+ 



+ 



tATn fl r X s-( V ) 



JO JO 



{u<X s -(v)Z n (s-)} 



Zn(s-) 



z + (l-z)X a -(y) z+(l-z)X s „(v)_ 
x M(ds,dz,du) 

tATn 

bX.-(v)- 1 [y(p)-'y(v)Z n (s-)]d8. 



Since X s ^ (v) >l/n for < s < T n , by a simple generalization of Theorem 2.1 
one can show the pathwise uniqueness holds for (4.7). Then, setting Zt = 
hmn-^oo Z n (t) we have 



(4.8) 



X t {p) = Z t X t (v)l {t<Too} + X t (v)l {t > Too} , t > 0. 
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Now from (4.7) and (4.8) we infer that {Xt(p) :t > 0} is measurable with 
respect to the cr-algebra & v generated by the process {X t (v) : t > 0} and the 
restricted martingale measures 

l {u<X s -{v)} B (ds,du), l {u < Xa _ {v)} M(ds,dz,du). 

By similar arguments, for any q S (y, 1] one can see {1 — Xt(q):t> 0} is mea- 
surable with respect to the cr-algebra % generated by the process {1 — X t (v)\ 
t > 0} and the restricted martingale measures 

l{x s -(v)<u<i}B(ds,du), l {Xs _ {v)<u < 1} M(ds,dz,du). 

Observe that {B(ds,X s _(v) + du)} is a white noise with intensity dsdu 
and {M(ds,dz,X s _(v) + du)} is a Poisson random measure with intensity 
dsv{dz)du. Then, given {X t (v) :t > 0} the cr-algebras & v and % are con- 
ditionally independent. That implies the Markov property of {(X(v),^ v ) : 
0<v<l}. □ 

Theorem 4.4. The path-valued Markov process {X(v) :0 < v < 1} has 
a p-cadlag modification. Consequently, there is a version of the solution flow 
{X t (v) : t > 0, < v < 1} of (4-3) with the following properties: 

(i) for each v £ [0,1], t^- X t (v) is cadlag on [0,oo) and solves (4-3); 

(ii) for each t>0, vt- > X t (v) is nondecreasing and cadlag on [0,1] with 
X t (0) > and X t (l) < 1. 

Proof. This follows from Lemma 3.5 and Proposition 4.2 by arguments 
as in the proof of Theorem 3.6. □ 

We call the solution flow {X t {v):t > 0,v £ [0,1]} of (4.3) specified in 
Theorem 4.4 a generalized Fleming -Viot flow following Bertoin and Le Gall 
(2003, 2005, 2006). The law of the flow is determined by the parameters 
(cr,6,7,z/). 

Let F[0, 1] be the set of nondecreasing cadlag functions / on [0, 1] such 
that < /(0) < /(l) < 1. Given a finite stopping time r and a function 
H G F[0, 1], let {X? t (v):t> 0} be the solution of 

X* t (v) = ii(v)+ [ T+t f 1 a[l {u < x » s _ {v)} -X? tS _(v)]B(ds,du) 

J T JO 

(4.9) + J T \ [l{v) - X ^ s _{v)]ds 

+ [ II 4Mu<x- s _(v)}-Ks-( v )W(ds,dz,du) 

J T JO Jo 

and write simply {X^{v) :t > 0} instead of {XQ t (v) :t > 0}. The pathwise 
uniqueness for the above equation follows from that of (4.3). Let F T) t be 
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the random operator on -F[0, 1] that maps /i to X^ t . As for the flow of 
CBI-processes we have 

Theorem 4.5. For any finite stopping time r we have P{X^ +t = F Tjt X^ 
for all t > 0} = 1. 

For any sub-probability measure fi(dv ) on [0, 1] with distribution func- 
tion v \— > [i(v), we write X^{dv) for the random sub-probability measure on 
[0, 1] determined by the random function v i— > X^(v). We call {Xj 1 : t > 0} the 
generalized Fleming-Viot process associated with the flow {X^{v) :t>0,v£ 
[0, 1]}. The reader may refer to Dawson (1993) and Ethier and Kurtz (1993) 
for the theory of classical Fleming-Viot processes. To give some character- 
izations of the generalized Fleming-Viot process, let us consider the step 
function 

n 

(4.10) f( U )=C l {0} (u)+J2 C i 1 (c H ^,a i ](u), ue[0,l], 

i=l 

where {co, ci, . . . , c n } C R and {0 = ciq < ai < ■ ■ ■ < a n = 1} is a partition of 
[0, 1]. For this function we have 

n 

(4.11) (XJt t f)=c o X?(0)+£c i [Xjt(a i )-X?(a i - 1 )]. 

i=i 

By (4.9) and (4.11) we have 

(X?, f) = (»J)+ f f a[g^(u)-(X^,f)}B(ds,du) 
Jo Jo 

(4-12) + fb^D-iX^J^ds 

Jo 

+ f f f z[g^{u) - {XZ_,f)]M(ds,dz,du), 
Jo Jo Jo 

where 

n 

(4.13) g${u) = c 1 {u < x m (0)} +J2 Cil {X?(<H-i)<u<XZ(a i )}- 

i=i 

The proofs of the following three results are similar to those for CBI-processes. 

Theorem 4.6. The generalized Fleming-Viot process {Xj* :t > 0} de- 
fined above is an almost surely cadlag strong Markov process with X^ = [i. 

Proposition 4.7. For any t > and f £ B[0, 1] we have 

(4.14) E[(ATJ)] = (n,f)e- bt + < 7 ,/>(l - e~ ht ). 
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Theorem 4.8. For any feB[0, 1] the process {(A"f , /) :t>0} has a cadlag 
modification. Moreover, there is a locally bounded function 1 1— > C(t) so that 

(4.15) e[su P (XZJ)} <C(t)[( M ,/} + ( 7) /} + {/i > / 2 } 1/2 + (7 J / 2 ) 1/2 ] 

for anyt>0 and f G B[0, 1] + . 

The generalized Fleming- Viot process can be characterized in terms of 
a martingale problem. Given any finite family {/i, . . . , / p } C -B[0, 1], write 



(4.16) 



G! P ,{/0( 7 ?) = fl( 7 ?./i)> r?GM![0,l]. 



i=l 



Let @i(L) be the linear span of the functions on Mi[0, 1] of the form (4.16), 
and let L be the linear operator on S'i(L) defined by 



i<j 



(4.17) 



/C{l,...,p},|/|>2 
V 



i=l 



fc=l 



fc=i 



(7,/i>n^/*)-n^/* 



fe=l 



where |/| denotes the cardinality of I C {1, . . . ,p} and 

Jo 

Theorem 4.9. TTie generalized Fleming-Viot process {Xj*:t>0} is 
the unique solution of the following martingale problem: for any p>l and 
{/i,-,/p}cB[0,l], 

(4.18) G pAfi} (X?) = G pm {n) + f LG P){ft} (X?) ds + mart. 

Jo 

Proof. We first consider a collection of step functions {/i, . . . , f p }. Let 
g^(s,u) be defined by (4.13) with / = /j. Since the compensation of the 
Poisson random measure in (4.12) can be disregarded, by Ito's formula we get 

G P,{fi}( X t) 



du 
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+ / ds [ v{dz) [ \f[[(X?J k )+zh$(8,u)]-f[(X?,f k )\du 
Jo Jo Jo [ k=1 k=l J 

ft p 

+ b / £[(7, fi) ~ (X»,fi)] IIW, /*) ds + mart. 

J° i=l k^i 

Gp^yM + o 2 fds /VlfHI^u) H (X?,f k )]x?(du) 

+ j* d S 1 jT 1 j j] [(x^, / fc ) + <(«)] - n ( x sJk) \x?(du) 



WE 



t=i 



(is + mart., 



(7,/i>iiw.M-nw^ 

fc^j fe=i 

where hf(s,u) =g?(s,u) - {Xg,fi) and Zf(tt) = fi(u) - (Xg,fi). It is simple 
to show 

1 l^{u)l^u)X^du) = (XjtJift) - (XjtJMXZJi). 
Then we continue with 



G P,{ft}( X t) 



G P,{/ i }(/ i ) 



a 



«<j L kj^i,j k=l 



ds 



ds [ u{dz) 
'o Vo 



i r p 



,fc=l 



fe=i 



WE 



fc=i 



(is + mart. 



i<j L fe^ij fe=l 

+ f ds t\(dz) t 1 \ e ^(i-^n^nw'/j) 

■ /o ^ ^ l/c{i,..,p} ieJ i*/ 



fe=i 
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r t v 

WE 



o 



i=l 



(7,fi)ll(X»,f k )-ll(XZ,f k 



k^i 



k=l 



ds + mart. 



G P ,{h}^) + ° 2 I E 

- 70 i<j 

JO JO JO Tr - r , 



k=l 



ds 



/C{l,...,p} 



,16/ j£J 



t V 



tl(x?,h) 

k=l 
ds + mart. 



>X»(du) 



Jq i=l L fc^j fe=l 

That gives (4.18) for step functions {/i,...,/ p }. For {/i,...,/ p } CB[0,1] 
one can show (4.18) by approximating the functions in the space L 2 (fi + 7) 
using bounded sequences of step functions. Since {X 1 ^ : t > 0} is a Markov 
process, and S>\{L) separates probability measures on M[0, 1], the unique- 
ness for the martingale problem holds [see Ethier and Kurtz (1986), page 182]. 
□ 

In particular, if /z(l) = 7(1) = 1, we have X^(l) = 1 for all t > 0, and the 
corresponding generalized Fleming- Viot process {X^ : t > 0} is a probability- 
valued Markov process with generator L defined by 

p 



LG p , {h] {r ] )=<J 2 Y; 



i<j 



(4.19) 



+ E 

JC{l,...,p},|/|>2 
V 

+E^Afi>n^>/fc 



k=l 

iei ' j0 k=l 



i=l 



where 



k^i 



Af(x) = b / [f(y)-f(x)h(dy), a: €[0,1]. 
'[0,1] 



This is a generalization of a classical Fleming- Viot process [see, e.g., Ethier 
and Kurtz (1993), page 351]. On the other hand, for 6 = the solution flow 
{Xt(v) :t > 0, < v < 1} of (4.3) corresponds to the A-coalescent process 
with A(dz) = a 2 5o + z 2 v(dz\ which is clear from (4.18) and the martingale 
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problem given by Theorem 1 in Bertoin and Le Gall (2005). For b > it seems 
the flow determines a coalescent process with a spatial structure. A serious 
exploration in the subject would be of interest to the understanding of the 
related dynamic systems. 

5. Scaling limit theorems. In this section, we prove some limit theorems 
for the generalized Fleming- Viot flows. We shall present the results in the 
setting of measure-valued processes and through the use of Markov process 
arguments. These are different from the approach of Bertoin and Le Gall 
(2006), who used the analysis of characteristics of semimartingales. For each 
k > 1 let cjfc > and b k > be two constants, let z 2 Vk(dz) be a finite mea- 
sure on (0,1] and let v t— > Jk( v ) be a nondecreasing continuous function 
on [0, 1] so that < 7fc(«) < 1 for all < v < 1. We denote by jk(dv) the 
sub-probability measure on [0, 1] so that 7fc([0, v]) = 7fc(v) for < v < 1. Let 
{Xt(v ) :t>0,v£ [0, 1]} be a generalized Fleming- Viot flow with parameters 
(o"fc,frfc,7fc,^fc) an d with X$(v) = v for v £ [0, 1]. Let Y k (t,v) = hX% t {k~~ l v) 
for t > and v £ [0, A;]. Let r]k(z) = kjk{k~ 1 z) and mk(dz) = Vk(k~ l dz) for 
z £ (0, k]. In view of (4.3), we can also define {Yk(t,v) :t>0,v£ [0, k]} di- 
rectly by 



where {Wk(ds, du)} is a white noise on (0, oo) x (0,k] with intensity dsdu, 
and {Nk(ds, dz, du)} is a Poisson random measure on (0,oo) x (0,k] 2 with 
intensity dsrrik{dz) du. In the sequel, we assume k>a for fixed a constant 
a > 0. Then the rescaled flow {Yfe(i, v) : t > 0, v £ [0, k]} induces an M[0, a\- 
valued process {Y£(t) :t > 0}. We are interested in the asymptotic behavior 
of {Y£(t) : t > 0} as k — > oo. Recall that A denotes the Lebesgue measure on 



(5.1) 




[0,oo). 
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1 r l 

1 '„2_2 / /-tll/iv-at 



+ I G"((Y h a (s),f))(Y k a (s),nd S 



where 



+ ds m k (dz) {G({Y k a (s),f) + zf(x))-G({Y k a (s)J)) 
JO JO J[0,a] 

-G'((Y k a ( S ),f))zf(x)}Y k a (s,dx) 

rt pk 

+ ds [ek{s,z) + £k{s,z)]mk{dz) + local mart, 
Jo Jo 



e k (s,z) = I k {G((Y k a ( s )J) + z[f(x) - k-\Y k «(s)J)]) 

-G((Y k a (s),f)+zf(x)) 
-k-'G'dY^JMY^J^Y^dx) 



and 



£k(s,z) = [k-Y k (s,a)\ 

x[G((Y k a ( S )J)-k- 1 z(Y k a (s)J)) 

- G((Y k a (s), /)) + k^G'((Y k a (s),f))z(Y k a (s), /)]. 

Proof. For the step function defined by (3.13) we get from (5.1) that 

(Y k a (t),f) = (X,f) + ka k f [ k h k (s-,u)W k (ds,du) 

Jo Jo 

(5-2) +kb k [\{vk,f)-{Y£(8-),f)]ds 

Jo 

pt pk pk 

+ / / / zh k (s-,u)N k (ds,dz,du), 



'o Jo Jo 

-1 Ivai 



where h k (s,u) =g k (s,u) - k (Y k (s),f) and 

n 

(5.3) 9k(s,U) = C l{u<Y k (s,0)} + ^2 c i 1 {Y k (s,a^ 1 )<u<Y k (s,a l )}- 

i=l 

Let l k (s,x) = f(x) - fc" 1 (*£(*),/). By (5.2) and Ito's formula, 
G((Y k a (t)J)) 

= G((X, /)) + kb k f G'((Y k a (s)J))[(rj k , f) - {Y k a (s),f)] ds 
Jo 
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+ 2***/ G"((Y k a (s),f))ds I h k (s,u) 2 du 



+ ds m k (dz) / {G((Y k a (s),f) + zh k (s,u)) 
Jo Jo Jo 

-G((Y k a (s),f)) 
-G'((Y k a (s),f))zh k (s,u)}du 

+ local mart. 

G«A, /)) + kb k f G'((Y k a (s)J))[( m , f) - (Y k a (s),f)} ds 



+ ^k 2 a 2 k I G"((Y k a (s),f))[(Y k a (s),n-k~t(Y k a (s),f}}ds 



1 l\ro,i 



+ / ds I m k {dz) I {G({Y k a (s),f) + zl k (s,x)) 

[0,o] 



Jo 

-G((Y k a (s)J)) 

-G'((Y k a (s),f))zl k (s,x)}Y k a (s,dx) 

+ [\k-Y k (s,a)]ds [ {G((Y k a (s),f)-k- 1 z(Y k a (s)J}) 
Jo Jo 

-G({Y k a (s),f)) 

+ k^G\{Y k a {s\ f))z(Y k a (s),f)}m k (dz) 

+ local mart. 

That gives the desired result for the step function. For / £ C[0, a] it follows 
by approximating the function by a sequence of step functions. □ 

Lemma 5.2. For t>0 and f £ C[0,a] + we have 



E 



sup (Y k a (s),f) 

0<s<t 



<(X,f) + kb k ( Vk ,f)t + U[(X,f) + (r ]k ,f)] I zm k {dz) 
+ 2v^[(A,/ 2 ) + (r ?fc ,/ 2 )] 1 /2 



1 \ 1/2' 

z 2 mk{dz 



o 



Proof. We first consider a nonnegative step function given by (3.13) 
with {cq, ci, . . . , c n } C R+. Let gt(s, u) and /i/t(s, u) be defined as in the proof 
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E 



sup (Y k a (s),f) 

0<s<t 



<(X,f) + 2ka k E 1 / 2 
+ kb k (r] k ,f)t + B 

»t rk rk 



V ft rk 



JO 



-i 2 



h k (s—,u)W(ds, du) 



o 



ds / zm k (dz) I \hk(s— , u)\ du 



o 



+ E 



Jl JO 



+ 2E 1/2 

It then follows that 

sup (Y k a (s),f) 

0<s<t 



z\h k (s—,u)\N k (ds, dz, du) 

t rl rk 

o Jo Jo 



zh k (s—,u)N k (ds, dz, du) 



< (X,f) + 2fco- fc E 1/2 |^ ds J h k (s,u) 2 du 



( rt rk rk ^ 

+ kb k {rj k J)t + 2~E<yJ ds J zm k (dz) J \h k (s, u)\ duj 



z 2 m k (dz) / h k (s,u) 2 du 



<(\J) + kb k (rj k J)t + 4E 



+ 2E 1/2 



(Y k a (s),f 2 )ds 



o 

(Y k a (s),f)ds j\m k {dz) 

1 \ 1/2 



kcr k + yj z 2 m k (dz] 



By Proposition 4.7 one can see 

n(Y k a (t), /)] = (A, f)e- kb ^ + ( m ,f)(l - e- kb ^) < (A, /) + ( Vk , /). 

Then we have the desired inequality for the step function. The inequality for 
/ G C[0, a] + follows by approximating this function with a bounded sequence 
of positive step functions. □ 

Lemma 5.3. Let r k be a bounded stopping time for {Y k (t) : t > 0}. Then 
for any t>0 and f £ C[0,a] we have 

B{\(Y k a (r k + t),f)-(Y k a (r k ),f)\} 



< E i/2 



(Y k a (r k +s),f 2 )ds 



1 \ 1/2 

ka k + ( / z 2 m k (dz) 
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f\( Vk ,\f\) + (Y k a (r k + s),\f\))ds 
Jo 



{Y k a (T k + s),\f\)ds / zm k (dz) 



Proof. We first consider the step function given by (3.13). Let g k (s,u) 
and h k (s,u) be denned as in the proof of Lemma 5.1. From (5.2) we have 

E{\(Y k a (r k + t),f)-(Y k a (r k )J)\} 
( r rt rk 

<ka k E 1 ' 2 < 



+ kb k E 

+ E l/2 

+ E 
+ E 



nh k (r k + s-,u)W(r k + ds, du) 

\(m,f)-(Y k a (r k + s-)J)\ds 
-t r l rk 

zh k (r k + s-, u)N k (r k + ds, dz, du) 




o Jo Jo 

t rk rk 



z\h k (r k + s-,u)\N k (r k + ds,dz,du) 



Jl Jo 

t rk rk 

ds I zm k (dz) / \h k (r k + s—,u)\du 
o Ji Jo 

By the property of independent increments of the white noise and the Pois- 
son random measure, 

E{\(Y k a (r k + t),f)-(Y k a (r k ),f)\} 



< A;cr fc E 1 / 2 |y ds J h k (r k + s,u) 2 du 



kb k E 



({Vk,\f\) + {Y k a (r k + s),\f\))ds 



+ E 1/2 { f ds f z 2 m k (dz) f h k {r k + s,u) 2 du 
[Jo Jo Jo 



+ 2E 



< E i/2 



ds / zm k (dz) / \h k (r k + s,u)\du 



(Y k a (r k + s),f 2 }ds 



ka k + 



1 \ V2' 

z 2 m k (dz) 



kb k ~E 
4E 



((Vk,\f\) + (Y k a (r k + s),\f\))ds 



{Y k a (T k + s),\f\)ds / zm k (dz) 
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Then (5.4) holds for the step function. For / £ C[0, a] the inequality follows 
by an approximation argument. □ 

Lemma 5.4. Suppose that kb k — > b, rj k — ?• n weakly on [0, a] and A; 2 a 2 , x 
<5o(dz) + (2 A z 2 )m k {dz) converges weakly on [0, oo) to a finite measure a 2 x 
5o(dz) + (2 A z 2 )m(dz) as k — > 00. Let {0 < ai < • • • < a n } 5e an ordered set of 
constants. Then {(Y^ 1 (t), . . . , Y? n (t)) : t > 0}, k = 1,2, . . . is a tight sequence 
in D([0,oo),M[0,oi] x ••• x M[0,a n ]). 

Proof. Let T k be a bounded stopping time for {Y£(t) :t > 0} and as- 
sume the sequence {r^ : k = 1,2, . . .} is uniformly bounded. Let fi S C[0, 
for i = 1, . . . , n. By (5.4) we see 

v\it\( Y k ai (n + t)Ji)- (Yt (T k ),fi)\) 



< 



/2 



(5.5) 



i=l 



{Y^{r k + s),f 2 )ds 



ka k + 



1 \ 1/2 

2 



z m k (dz) 



n [ ft 
i=i L- 70 



+ 4 E E 



i=l 



( 1 fc i ( r fc + s )»l/il) rf ' s / zm k {dz 



Then the inequality in Lemma 5.2 implies 



limsupE^|(yf (n + t)Ji) - (^(r fe ),/i)| =0. 

By a criterion of Aldous (1978), the sequence {((Y fc ai (t),/i), . . . , jfn)) : 

t > 0} is tight in D([0, 00), M n ) [see also Ethier and Kurtz (1986), pages 
137 and 138]. Then a simple extension of the tightness criterion of Roelly- 
Coppoletta (1986) implies {{Y* 1 (t), Y£" (t)) : t > 0} is tight in D([0, 00), 
M[0,ai] x ■•■ x M[0,a n ]). □ 



Suppose that a > and 6 > are two constants, u 1— > n(-u) is a nonnega- 
tive and nondecreasing continuous function on [0,oo) and (z A z 2 )m(dz) is 
a finite measure on (0,oo). Let n(dv) be the Radon measure on [0, 00) so 
that rj([0, v]) = i](v) for v > 0. Suppose that {W(ds,du)} is a white noise on 
(0,oo) 2 with intensity dsdz and {N(ds, dz, du)} is a Poisson random mea- 
sure on (0,oo) 3 with intensity dsm(dz) du. Let {X t (v) :t > 0,v > 0} be the 
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solution flow of the stochastic equation 

ft rX s -(v) rt 

X t (v) = v + a / W(ds,du) + b [n(v)-X s _(v)]ds 



o Jo 



(5.6) 

rt poo rX a —(v) 



lo Jo Jo 



zN(ds, dz, du). 



By Theorem 3.11, for each a > the flow {Xt(v) :t > 0,v > 0} induces an 
M[0, a] -valued immigration superprocess {Xf :t > 0} which is the unique 
solution of the following martingale problem: for every G G C 2 (]R) and / € 
C[0,a], 

C((X t ,/}) 

= G((X,f))+b [ G / ((X s ,/))[(r ? ,/)-(X s ,/)]d S 
■/ o 

+ l<T 2 J*G"((X s ,f))(X s ,f)ds 



(5.7) 



+ / ds / m(dz) / [G((X s J) + zf(x)) 
Jo Jo 

- G((X S J)) - G'((X s ,f))zf(x)]X s (dx) 

+ local mart. 



Theorem 5.5. Suppose that kb^ — > b, rjk rj weakly on [0,a] and 
k 2 a?6o(dz) + (z A z 2 )mk(dz) converges weakly on [0, oo) to a finite measure 
a 2 5o(dz) + (z A z 2 )m(dz) as k — >■ oo. Then {Y^(t) :t > 0} converges to the 
immigration superprocess {Xf 1 :t > 0} in distribution on D([0,oo),M[0,a]). 

For the proof of the above theorem, let us make some preparations. 
Since the solution of the martingale problem (5.7) is unique, it suffices to 
prove any weak limit point {Z®:t > 0} of the sequence \Y£(t):t > 0} is 
the solution of the martingale problem. To simplify the notation we pass to 
a subsequence and simply assume {Y£(t):t > 0} converges to {Zf :t > 0} 
in distribution. Using Skorokhod's representation theorem, we can also as- 
sume \Y£{t):t > 0} and {Zf:t > 0} are defined on the same probability 
space and {Y£(t):t > 0} converges a.s. to {Z^:t > 0} in the topology of 
D([0,oo),M[0,a]). For n> 1 let 

r„ = inf(t>0:sup / [1 + (Y£(s) + Z s a , l) 2 ] ds > n\ . 

[ k>l Jo J 

It is easy to see that r n — > oo as n — > oo. 
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Lemma 5.6. Suppose that kb k —>b, rj k V weakly on [0, a] and k 2 o~\ x 
6o(dz) + (z A z 2 )nik{dz) converges weakly on [0, oo) to a finite measure a 2 x 
So(dz) + (z A z 2 )m(dz) as k — > oo. Let e k (s, z) 6e defined as in Lemma 5.1. 
Then for each n>\ we have 



E 



tAr„ 



ds / \e k (s,z)\m k (dz) 



0. 



k oo. 



Proof. By the mean- value theorem, we have 
e k (s,z) = ±z(Y k a (s)J) 



x / [G'({Y k a (s),f) + z0 k (s,x)) -G'((Y k a (s),f))]Y k a (s,dx), 
Jo 

where 9 k (s,x) takes values between f(x) and f(x) — k~ 1 (Y k (s), /). Conse- 
quently, 



z(n a ( S ),i) 2 . 



|e*(*,«)l<f||^IW(«).l/l>WW,l><f n/ 
Moreover, since (Y k (s), 1) < fe, we get 

\e k (s,z)\<l\\G"\\z 2 (Y k a (s),\f\) f k \e k (s,x)\Y k a (s,dx) 



<-||G"||^(y fc a ( S ),|/|) [\f(x)\+k^(Y k a (s),\f\)}Y k a (s,dx) 



<-imi 2 ||G"||z 2 <Y fc %<0,l> 2 



It follows that 



E 



tAT„ 



tAr n 



(Y k a (s),l) 2 ds 



ds / \e k {s,z)\m k (dz) 
Jo 

C f k 

< — I (z A z 2 )m k (dz)E 
k Jo 

where C = 2||/||(||G'|| + ||G"||||/||). The right-hand side goes to zero as k 
oo. □ 



Lemma 5.7. Suppose that kb k — > b, rjk—tr] weakly on [0, a] and k 2 a k x 
5o(dz) + {z A z 2 )m k (dz) converges weakly on [0, oo) to a finite measure a 2 x 
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So(dz) + (z A z 2 )m(dz) as k—>oo. Let £ k (s,z) be defined as in Lemma 5.1. 
Then for each n > 1 we have 



E 



tAr n 



UO 



ds / \£ k (s,z)\m k (dz) 



0. 



k—^oo. 



Proof. It is elementary to see that 

Ms, z)\ < k\G({Y h a (s), f) - fc-^(y fc a ( a ), /)) 

- G((Y k a (s), /)) + k~ l G'((Y k a (s),f))z(Y k a (s), f)\ 

< mm^WG'MY^s), |/|), ±-\\G"\\z 2 (Y k «(s), |/|) 2 

<C[l + (Y fc >),l> 2 ](*A*rV), 

where C= ||/||(2||G'|| + ||/||||G"||/2). Then we have 

rt/\T n 

e' 



ds I \£k{s,z)\m k {dz 
Jo 

rk ( rthTn, 

<C J (zAk- 1 z 2 )m k (dz)Elj [l + (Y k a (s),l) 2 ]ds 



<nC (z Ak~ l z 2 )m k {dz). 
Jo 

The right-hand side tends to zero as k — > oo. □ 

PROOF of Theorem 5.5. Let / e C[0,a]. Then {(Y k a (t), f) :t > 0} con- 
verges a.s. to {(Zf,f) :t > 0} in the topology of D([0,oo),R). Consequently 
we have a.s. (Y k (t),f) — > (Z£,f) for a.e. t > [see, e.g., Ethier and Kurtz 
(1986), page 118]. By Lemma 5.1, 

G((Y k a (t), /)) = G«A, /)) + fc6 fe / G"«Y fe a ( s ), /))(%,/) 



/ G'((F fc a ( S ),/))(Y fc a ( S ),/)d S 

1 



+ ^k 2 a 2 I G"((Y k a (s),f))(Y k a (s),nds 



(5.8) 



~\hJ q G"((Y k a (s),f))(Y£(s)J) 2 ds 

+ f ds f k m k (dz) [ H(x,z,(Z*J))Y k a (s,dx) 

JO Jo J\0,a] 
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where 



and 



+ ds [e k (s,z) + £ k (s,z) + Ck(s,z)]m k (dz) 
Jo Jo 

+ local mart., 



H(x, z, u) = G(u + zf(x)) - G(u) - G'{u)zf{x) 



Ck(s,z)= / [H(x,z,(Y k a (s),f))-H(x,z,(Z^f))]Y k a (s,dx). 

J [0,o] 

By the mean- value theorem, 

\C k (s,z)\< [ \H' u (x,z,6 k (s))(Y k a (s) - Z«,f)\Y k a (s,dx), 

J[0,k] 

where 9 k (s) takes values between (Y k (s),f) and (Z°,/). For G G C 3 (R) we 
have 



\H' u (x,z,9 k (s))\ = \G'(9 k (s) + zf(x)) - G'(9 k (s)) - G"(9 k (s))zf(x)\ 
<||/||(2||G"|| + i||/||||G'"||)(,Az 2 ). 

It follows that 

|a( S )|<||/||(2||G"|| + i||/||||G w ||)(.Az 2 ) 
x(Y k a (s),l)\(Y k a (s)-Z«,f)\. 
By (5.9) and Schwarz's inequality, 



(5.9) 



E 



tAT n 



UO 



ds / \( k (s)\m k (dz) 



< C k (t) \ E 



t/\T„ 



(Y k a (s)-Z^f) 2 ds 



1/2 



x <^ E 



t/\T n 



(Y k a (s),lfds 



<V^C k (t)\E 



tAT n 



1/2 



1/2 



where 



C7 fc (t) = ||/||(2||Gr|| + i||G w 



[z A z 2 )m k (dz). 
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Note that sup fc>1 Ck(t) < oo. It then follows that 



E 



tAT n rk 

ds / \Ck(s)\m k (dz) 



JO 



0, k— >oo. 



Now letting k — > oo in (5.8) and using Lemmas 5.6 and 5.7 we obtain (5.7) for 
GgC 3 (R) . A simple approximation shows the martingale problem actually 
holds for any Ge C 2 {R). □ 

Theorem 5.8. Suppose that kb k — > b, n k —> n weakly on [0, a] and 
k 2 a 2 5o(dz) + (z A z 2 )mk(dz) converges weakly on [0,oo) to a finite mea- 
sure a 2 5o(dz) + (z A z 2 )m(dz) as k — >■ oo. Lei {0 < a\ < ■ ■ ■ < a n = a} 6e 
an ordered set of constants. Then {(Y^ 1 (t), . . . , Y" fc a ™ (t)) :t > 0} converges to 
{(X? 1 , . . . ,X t a ") :t > 0} m distribution on D([0, oo), M [0, Oi] x • • • x M[0, a n ]). 

Proof. By Lemma 5.4 the sequence {(y fe ai (*),..., Y" fc a " (t)) : t > 0} is tight 
in £>([0,oo), M[0,ai] x • • • x M[0,a n ]). Let {(Zf 1 ,.- • :t > °l be a weak 

limit point of {(Y^(t), . . .,Y£ n (t)):t> 0}. To get the result, we only need 
to show {(Z? 1 , . . . , Zf") : t > 0} and {(Xf 1 , . . . ,X^) :t > 0} have identical 
distributions on D([0, oo), M[0, ai] x • • • x M[0,o n ]). By passing to a subse- 
quence and using Skorokhod's representation, we can assume {(Y" fc ai (t), . . . , 
Y£ n (t)):t > 0} converges to {(Z t ai , . . . , Zf n ) : f > 0} almost surely in the 
topology of D([0,oo), M[0,oi] x • • • x M[0,a n \). Theorem 5.5 implies {Z? n : 
t > 0} is an immigration superprocess solving the martingale problem (5.7) 
with a = a n . Let Z^ denote the restriction of Z^ n to [0, a*]. Then Z? n = Z? n 
in particular. We will show {(Zf 1 , . . . , Zf") : t > 0} and {(Zf 1 , . . . , Z? n ) : t > 0} 
are indistinguishable. That will imply the desired result since {(X^ 1 , . . . , 
X? n ) :t>0} and {(Z^ 1 , . . . , Z^ n ) : t > 0} clearly have identical distributions 
on -D([0, oo), M[0, a\] x • • • x M[0, a n ]). By the general theory of cadlag pro- 
cesses, the complement in [0, oo) of 

D(Z) -={t>0: P(Z t fll = Z?± , . . . , Z? n = Z£) = 1} 

is at most countable [see Ethier and Kurtz (1986), page 131]. For any t G 
D(Z) we have almost surely lim^oo Y£' (t) = Zp for each i = 1, . . . , n [see 
Ethier and Kurtz (1986), page 118]. By an elementary property of weak 
convergence, for any t G D(Z) we almost surely have 

Zp (M)= lim (t,[0,Oi])= lim Y^(t,[0, ai ]) 

k— >oo re— >-oo 

< 2^([0,Oi]) = ^([O.Oi]) = Zf([0,a,]). 
Since Theorem 5.5 implies {Zp : t > 0} is equivalent to {Zp : t > 0}, we have 
E[Zf([0,a i ])]=E[Zf([0 ) o i ])]. 
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It then follows that almost surely 

(5.10) lim Y?(t,[0,a i ])=Z?([0,a i ]). 

On the other hand, since Y£ n {t) — > Z^ n , for any closed set C C [0, Oj] we 
have 

(5.11) lim sup (t,C) = lim Y* n (t,C) < Z? n {C) = Z?(C). 

k— >oo fe— >oo 

By (5.10) and (5.11) we have = lim*^ (t) = Z^ . Then {Z? : t > 0} 
and {Zp : t > 0} are indistinguishable since both processes are cadlag. □ 

Let ^# be the space of Radon measures on [0, oo) furnished with a metric 
compatible with the vague convergence. The result of Theorem 5.8 clearly 
implies the convergence of {Yk(t) :t > 0} in distribution on D([0,oo),^#) 
with the Skorokhod topology. From Theorem 5.8 we can also derive the 
following generalization of a result of Bertoin and Le Gall (2006) [see also 
Bertoin and Le Gall (2000) for an earlier result]. 

Corollary 5.9. Suppose that kb^ — > b, weakly on [0,a] and 

k 2 o~?So(dz) + (z A z 2 )mk(dz) converges weakly on [0, oo) to a finite mea- 
sure a 2 5o(dz) |(zA z 2 )m{dz) as k — > oo. Let {0 < a\ < • • • < a n } be an 
ordered set of constants. Then {(Yk(t, a\), . . . , Yk(t, a n )) : t > 0} converges to 
{(Xt(ai),...,Xt(a n )):t>0} in distribution on D([0, oo),R™ ). 
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